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Abstract 

Using two sets of high-precision Monte Carlo data for the two- 
dimensional XY model in the Villain formulation on square L x L lat- 
tices, the scaling behavior of the susceptibility x and correlation length 
£ at the Kosterlitz-Thouless phase transition is analyzed with empha- 
sis on multiplicative logarithmic corrections (lnL) -2r in the finite-size 
scaling region and (ln£) -2r in the high-temperature phase near crit- 
icality, respectively. By analyzing the susceptibility at criticality on 
lattices of size up to 512 2 we obtain r = —0.0270(10), in agreement 
with recent work of Kenna and Irving on the the finite-size scaling of 
Lee- Yang zeros in the cosine formulation of the XY model. By study- 
ing susceptibilities and correlation lengths up to £ ~ 140 in the high- 
temperature phase, however, we arrive at quite a different estimate of 
r = 0.0560(17), which is in good agreement with recent analyses of 
thermodynamic Monte Carlo data and high-temperature series expan- 
sions of the cosine formulation. 



PACS numbers: 05.50.+q, 75.10.Hk, 64.60. Cn, 11.15.Ha 



1 Introduction 



Ever since the seminal work of Kosterlitz and Thouless (KT) in 1973, 1,2 the 
two-dimensional (2D) XY model has been the subject of extensive experi- 
mental, analytical and numerical investigations. 3 Physically the interest in 
this model arises from studies of layers of superconducting materials and 
films of liquid helium, 4 Josephson-j unction arrays, 5 and some magnetic sys- 
tems. 6 Theoretically the peculiar behavior of the KT phase transition, which 
is believed to be driven by the unbinding of defect pairs, has attracted much 
interest. Despite all these efforts, however, the details of the phase transition 
are not yet fully understood. 

In a recent Monte Carlo (MC) simulation study of Lee- Yang partition 
function zeros, Kenna and Irving 7 ' 8 raised again the question of logarithmic 
corrections 2 ' 9 to the leading finite-size scaling (FSS) scaling behavior. If the 
linear lattice size is denoted by L and the multiplicative logarithmic correc- 
tions are parametrized as (lnL)~ 2r , their numerical result is r = —0.02(1), 
while the standard KT theory would predict quite a different exponent of 
r = —1/16 = — 0.0625. 2 ' 9 Moreover, by reanalyzing "thermodynamic" MC 
data of Refs. 10,11 obtained on lattices with L > 7£, where £ is the correlation 
length, Patrascioiu and Seiler 12 obtained an estimate of r = 0.077(46), and 
by analyzing long high-temperature series expansions, Campostrini et al. 13 
also arrived at positive values in the range r = 0.042(5) to r = 0.05(2), de- 
pending on the quantity considered. While the estimates of the latter two 
groups are consistent with each other, they are incompatible with the FSS 
result of Kenna and Irving, which, on the other hand, is somewhat "closer" 
to the theoretical prediction. 

All numerical estimates quoted above were obtained in the cosine formu- 
lation of the XY model. The purpose of this note is to add further evidence in 
one or the other direction by analyzing the logarithmic corrections in the Vil- 
lain formulation 14 of the XY model, which is actually (sometimes implicitly) 
the starting point of most if not all theoretical investigations. 
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2 Scaling predictions 



In the Villain XY model 14 the Boltzmann factor of the cosine formulation, 
B cos = Il;E,iexp [/3 cos cos(Vi0(cc))], is replaced by the periodic Gaussian 



B = H £ exp 

X.i n=— oo 



■|(V i fl-27m) 5 



(1) 



where /3 is the inverse temperature in natural units, and Vj# = 9(x + i) — 
9(x) are lattice gradients. A discussion of the relation between the two 
formulations as well as numerical comparisons can be found in Refs.. 14 ' 15 
The two-point correlation function (s— (cos(0), sin(0)) 



G(x) = (s(x) • s(0)) = (cos(9(x) - 9(0))) 
is predicted to behave at the critical temperature T c = 1/ f3 c as 9 

(\n\x\)- 2r 



G(x) oc 



l + O 



In In 


x\ 


In 


X 





(2) 



(3) 



with r = —1/16 and i] = 1/4. For the power of the logarithmic term we have 
adopted the notation of Refs.. 7,8 In the high-temperature phase near criti- 
cality, i.e. < t = T/T c -1<1, this implies for the magnetic susceptibility, 



a scaling behavior 



where 



X = V((Y,s(x)/V) 2 ) = Y,G(x), 

X X 



X oc£ 2 ^(lnO- 2 11 + 0(lnln£/lnO], 



(4) 



(5) 



f oc exp(br v ) (6) 

is the correlation length, with v — 1/2 and b being a non-universal positive 
constant. Expressing £ in terms of t, eq. (5) can also be written as 

X oc£ 2 -^[(1 + O(rint)]. (7) 

Very close to T c eq. (5) cannot hold for a finite system with linear size L <^ £. 
Here £ has to be replaced by L, and we expect to observe a FSS behavior 



X oc L 2 ^(lnL)- 2r [l + 0(mmL/mL)]. 



(8) 
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In numerical simulations it proved to be very difficult to verify the KT 
scaling laws unambiguously. However, if one rejects power-law ansatze with 
unnaturally large exponents and large confluent correction terms, then, among 
the two alternatives, a pure power-law or the exponential KT divergences, the 
KT predictions are clearly favored. This is the conclusion of most numerical 
studies of the cosine formulation 10 ' 11 and, with even stronger evidence, also 
of the Villain formulation 16 considered here. In this note we shall therefore 
not study this fundamental question again. We rather assume eqs. (5)- (8) 
to be qualitatively valid and try to determine the exponents 77, v, and r. 
Unfortunately even this goal is far too ambitious, since a precise determina- 
tion of all three critical exponents together with the (non-universal) value of 
f3 c would require much more accurate data than one can hope to generate 
with present day techniques. We therefore hold the exponents v — 1/2 and 
7] — 1/4 fixed at their theoretically predicted values and ask if any devia- 
tion of the data from the leading scaling behavior can be explained by the 
logarithmic corrections in eqs. (5) and (8). 

3 Results 

In Ref. 16 we have reported high-precision MC simulations of the Villain model 
(1), using the single-cluster update algorithm and improved estimators for 
the two-point correlation function. This enabled us to obtain on a 1200 2 
square lattice data for the correlation length up to £ ps 140. Since L > 8£ 
this value should be a very good approximation of the thermodynamic limit. 
By performing fits of £ to the KT prediction (6) and of x to (5) (without 
the logarithmic term) with four free parameters (the prefactor, b, u, and f3 c ) 
we obtained f3 c = 0.752(5) and v = 0.48(10). The estimate of f3 c is in very 
good agreement with the more precise value of f3 c = 0.7524(7) obtained in 
Ref. 17 from a study of the dual discrete Gaussian model (see also 18 ). Using 
the ansatz (7), i.e. including the theoretically predicted correction t -1 / 16 , did 
not improve the quality of the fits. 

Further data of the susceptibility at criticality on lattices with up to 
512 2 sites showed a clear scaling behavior for L > 100, x L 2 ~ v , with 
r] = 0.2495 « 1/4 at (3 = 0.74, and rj = 0.2389(6) ^ 1/4 at (5 = 0.75. This 
is obviously not consistent with the prediction that rj = 1/4 at f3 c . Since the 
estimate of f3 c from two completely independent simulations agreed so well we 
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concluded in Ref. 16 that i]((3 c ) ^ 1/4, in disagreement with the KT prediction. 
To reconcile simulations and theory we speculated that the scaling curve for 
X might still change for much larger system sizes, but this is of course not 
very convincing. Mainly based on our negative experience with the t^ 1 / 16 
correction in the x{T) fits, we did not try, however, to attribute the observed 
discrepancy to logarithmic corrections. 

The data at f3 — 0.75 and a fit in the range L > 64 according to 
ln(x/L 7 / 4 ) = const. + (1/4 — rj) InL is reproduced in Fig. 1(a). In Fig. 1(b) 
we show the same data, but now fix rj = 1/4 at the theoretical value 
and assume that (8) with the logarithmic correction is valid. Since then 
ln(x/L 7 / 4 ) = const. — 2rln(lnL), we expect a straight line when ln(x/L 7 / 4 ) 
is plotted against ln(lnL). As is demonstrated in Fig. 1(b) this is clearly the 
case. Also shown is a linear fit which is of high statistical quality (goodness- 
of-fit parameter Q = 0.61) and yields a slope of 0.0540(19), or 

r = -0.0270 ±0.0010, (9) 

in good agreement with the estimate of r = —0.02(1) from the FSS of Lee- 
Yang zeros in Refs.. 7 ' 8 To summarize this subsection, by allowing for logarith- 
mic corrections we can reconcile the numerical estimate of f3 c ~ 0.752 with the 
KT prediction i]((3 c ) = 1/4. The value of the exponent r, however, is clearly 
not in agreement with the theoretical prediction r = —1/16 = —0.0625. 

Let us next consider the scaling behavior of the thermodynamic data near 
criticality in the high-temperature phase. In Ref. 16 we neglected logarithmic 
corrections in (5) and tested the relation lnx/^ 7 ^ 4 = const. + (1/4 — 77) In £ 
graphically. This plot is reproduced in Fig. 2(a). We see that the curve has 
a negative slope, corresponding to 77 > 1/4. We also observe, however, that 
the data are curved and that for large £ the slope decreases. By defining i] cS 
from the local slopes, we obtained at the scale of £ ^ 110 . . . 140 an estimate 
of f] cS fx 0.267. Notice that this effective rj is above 1/4, while from FSS 
without logarithmic corrections we would have extracted an effective 77 that 
is smaller than 1/4. In Fig. 2(b) we show the same data, but similar to Fig. 1 
we now again fix 77 = 1/4 at the theoretical value and assume that (5) with the 
logarithmic correction is valid. Since then ln(x/£ 7 ^ 4 ) = const. — 2rm(ln£), 
we expect a straight line when ln(x/£ 7 ^ 4 ) is plotted against ln(m£). This is 
indeed the case, and from the fit over all available data points (with Q = 0.97) 
we obtain 

r = 0.0560 ±0.0017, (10) 
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in qualitative agreement with the results in Refs., 12,13 which are also derived 
from the approach to criticality in the high-temperature phase. The value 
(10) is clearly different from (9), and is very far from the theoretical estimate 
r = —1/16 = —0.0625. In retrospective this "explains" why we did not 
observe any improvement when trying fits of x{T) with the t r correction 
fixed to the theoretical prediction t -1 / 16 . 

We repeated the analysis leading to the Villain model estimate (10) also 
with the three data points for the cosine model in Ref. 16 (with £ ps 21, 40, 
and 70) and obtained a compatible value of r = 0.047(8). Furthermore, using 
the more extensive data sets of Ref. 10 and Ref. 11 we find compatible values 
of r = 0.050(10) and r = 0.049(10), respectively. 

We also tried to use the scaling form (7) which requires as input infor- 
mation the value of f3 c . Using the most accurate estimate of (3 C = 0.7524 we 
find the result shown in Fig. 3. Again the linear scaling looks almost perfect, 
but from the slope we now read off an even larger value of r = 0.0922(28). 
Qualitatively this can be understood as follows. Going from (5) to (7) we 
replace ln£ by t~ v = t~ l l 2 . Asymptotically this follows from the scaling be- 
havior of £ in (6). This implicitly assumes, however, that the constant in 
the proper relation, ln£ = const. + bt~ u , can be neglected. If t is not really 
asymptotically small, this is not justified. In fact, the plot of ln(ln£) vs — lnt 
in Fig. 4 does show effectively an almost linear behavior, but with a slope 
completely different from the asymptotic value v — 1/2. 

Finally it was of course tempting to enquire if the observed discrepancies 
between the numerical data and the theoretical expectations can be blamed 
on the additive logarithmic corrections in (5) and (8). To test this possibility 
we have replotted in Fig. 5 the data at criticality in the form x/L 2_T '(lnL)~ 2r 
vs In In L/ In L and the thermodynamic data in the form x/£ 2 ~ T '(ln£)~ 2r vs 
lnln£/ln£, assuming the theoretically predicted values of rj and r. The 
double valuedness in Fig. 5(b) is caused by the fact that f(x) = In In xj In x 
assumes a maximum / max = 1/e ~ 0.3679 at x max = e 6 ~ 15.15. We see 
that both the data for L > 64 or £ > 40 can be well fitted with a simple 
linear function. With a parabolic ansatz the acceptable fit range can even 
be extended to smaller values of L or £. From Fig. 5 it is obvious, however, 
that we are still too far away from the truly asymptotic region x — > to 
take this as a convincing evidence that additive logarithmic corrections can 
reconcile simulations and theory. 
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4 Discussion 



In summary we have shown that, when multiplicative logarithmic corrections 
are taken into account, numerical simulation data of the 2D XY Villain model 
are quite consistent with the leading KT predictions even at a quantitative 
level with critical exponents fixed to the theoretical values of v — 1/2 and 
i] = 1/4. Estimates of the logarithmic correction exponent r, however, turn 
out to be quite inconsistent. Scaling analyses in the FSS region yield a neg- 
ative (r ps —0.02 ... — 0.03) and analyses in the high-temparature phase a 
positive (r ps 0.04 . . . 0.08) value, both being quite different from the theoret- 
ical prediction of r = — 1/16 = —0.0625. This is obviously related to the fact 
that analyses neglecting logarithmic corrections tended to estimate rj > 1/4 
using thermodynamic data and rj < 1/4 in the FSS region. We have no good 
explanation for this observation other than the common, but unfortunately 
probably correct statement 19 ' 20 that the studied system sizes are still much 
too small to resolve these discrepancies. 
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Figure 1: Finite-size scaling of the susceptibility at criticality. If logarithmic 
corrections are neglected, the slope in (a) gives 1/4— 77. If 77 = 1/4 is assumed, 
the slope in (b) yields — 2r, the exponent of the logarithmic correction. 
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2D XY Villain 




ln(ln£) 



Figure 2: Test of the scaling relation x °c £ 2 T '(ln£) 2r in the range £ as 
10 . . . 140, rewritten as ln(x/£ 7/4 ) = const. + (1/4 - 77) ln£ - 2rm(ln£)- The 
linear behavior in (b) shows that the data are compatible with 77 = 1/4. 
As is already obvious from (a), the exponent r must then be positive, in 
disagreement with the theoretical prediction r = —1/16. 
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Figure 3: Test of scaling similar to Fig. 2, but with ln£ replaced by t = 
T/T c - 1 (cp. eq. (6)). 
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Figure 4: Correlation length vs reduced temperature. In the range £ ~ 
10 . . . 140 the slope is effectively about 0.8, while asymptotically it should 
approach v = 0.5 according to eq. (6). 
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Figure 5: Test for additive logarithmic corrections in (a) the data at criticality 
and (b) the thermodynamic data. Here the exponents rj and r are assumed 
to take the theoretically predicted values rj — 1/4 and r = —1/16. 
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